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Escaping from Modified Newtonian Dynamics 
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ABSTRACT 

We present a new test of modified Newtonian dynamics (MOND) on galactic scales, 
based on the escape speed in the Solar Neighbourhood. This test is independent from 
other empirical successes of MOND at reproducing the phenomenology of galactic 
rotation curves. The galactic escape speed in MOND is entirely determined by the 
baryonic content of the Galaxy and the external field in which it is embedded. We 
estimate that the external field in which the Milky Way must be embedded to produce 
the observed local escape speed of ~ 550 kms -1 is of the order of an/100, where ao 
is the dividing acceleration scale below which gravity is boosted in MOND. This is 
compatible with the external gravitational field actually acting on the Milky Way. 
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1 INTRODUCTION 

The speed needed by a star in the Solar Neighbourghood 
to escape the local gravitational field of the Milky Way - 
known as the local galactic escape speed - has long been 
known to be an essential indicator of the amount of mass 
beyond the galactocentric radius of the Sun (e.g. Binney 
Tremaine 1987). Given the poor constraints on the shape of 
the galactic rotation curve beyond that radius (e.g. Binney 
& Dehnen 1997; Famaey & Binney 2005, herefater FB05), it 
has become one of the most convincing arguments in favour 
of a halo of dark matter in our Galaxy (Smith et al. 2007). 
However, the observed tight correlation between the mass 
profiles of baryonic matter and dark matter at all radii in 
spiral galaxies (e.g. McGaugh et al. 2007; Famaey et al. 
2007) rather lends support to modified Newtonian dynam- 
ics (MOND). This theory postulates that for accelerations 
below ao = 1.2 x 10 -10 ms -2 the effective gravitational at- 
traction approaches (gNdo) 1 ^ 2 where grjv is the usual Newto- 
nian gravitational field. Without resorting to galactic dark 
matter, this simple prescription is amazingly successful at 
reproducing many aspects of the dynamics of galaxies that 
were unknown when Milgrom (1983) introduced the theory 
(e.g. Sanders & McGaugh 2002; Tiret & Combes 2007). It 
is however worth noting that the theory might still welcome 
(and might even need) hot dark matter, which could con- 
dense in galaxy clusters and boost the (modified) gravita- 
tional potential on Mpc scale (Sanders 2003; Angus et al. 
2007). 
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In MOND, the gravitational potential $ of an isolated 
point mass M at a distance r S> (GM/ao) 1 ^ 2 is of the form: 



k *(r) ~ (GMa ) 1/2 ln(r). 



(1) 



This implies that the associated rotation curve is asymptot- 
ically flat, but in turn it implies that there is no escape from 
such an isolated point mass. In practice however, no objects 
are truly isolated in the Universe and this has wider and 
more subtle implications in MOND than in Newton-Einstein 
gravity, for the very reason that MOND breaks the Strong 
Equivalence Principle (SEP). In particular, this letter illus- 
trates how the external field regularizes the above divergent 
potential, so that it is possible to escape from non-isolated 
point masses in MOND. 

More generally, this breakdown of the SEP remarkably 
allows to derive properties of the gravitational field in which 
a system is embedded from its internal dynamics (and not 
only from tides). For this reason, we are able to compute 
hereafter the external gravitational field in which the Milky 
Way should be embedded in MOND in order to account for 
the measured local galactic escape speed. The result is in 
good agreement with the gravitational field that Large Scale 
Structures are expected to exert on the Galaxy (of the order 
of oo/100). The MOND paradigm thus successfully explains 
the value of the local escape speed without referring to any 
halo of dark matter. 
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2 OBSERVED GALACTIC ESCAPE SPEED 

The local galactic escape speed in the Solar Neighbourhood 
can be estimated by using local samples of high velocity stars 
(Fich & Tremaine 1991), but this task was until recently 
hampered by the very small number of such observed stars 
and by the fact that they were selected from proper-motion 
surveys. Now, radial velocity data from the Radial Velocity 
Experiment (RAVE), a magnitude- limited survey avoiding 
kinematical biases (Steinmetz et al. 2006), have been used 
to constrain the local escape speed (Smith et al. 2007). A 
total of 33 stars (16 from the RAVE survey and 17 from 
Beers et al. 2000) with velocities higher than 300 kms -1 
have been selected, and, following the technique of Leonard 
& Tremaine (1990), the local escape speed has been esti- 
mated to lie in the following 90% confidence interval: 

498 km s -1 < v csc < 608 kms -1 , (2) 

with a median likelihood of 544 km s -1 . Any unbound star in 
the sample could have biased the result but Yu & Tremaine 
(2003) have shown that hypervelocity stars ejected from the 
galactic center are actually very rare (i.e. 1 every 10 5 years) 
and thus not likely to belong to such a survey. 



3 GALACTIC ESCAPE SPEED IN MOND 
3.1 The external field effect 

The local galactic escape speed essentially depends on the 
profile of the gravitational potential in the outskirts of the 
Milky Way. The latter may of course depend on the precise 
realization of the MOND phenomenology. Here we assume 
that the gravitational potential $ of the Milky Way satis- 
fies the modified Poisson equation (Bekenstein & Milgrom 
1984): 

V.KI - V$ + g c |/a )V$] ~4ttG*p, (3) 

where $ is the non-external part of the potential, p is the 
density of the (baryonic) matter, g e is the approximatively 
constant external field in which the Milky Way is embedded, 
and fJ.(x) is a function which runs from p(x) = 1 at x 3> 1 
to l-t(x) = x at x <C 1, and which can thus never be smaller 
than |g e |/ao. Note however that Eq. Q is only an approx- 
imate and effective way to take into account the effect of 
the external field on local physics, in order to avoid solv- 
ing the modified Poisson equation of Bekenstein & Milgrom 
(1984) with an external source term p e xt on the right-hand 
side. This type of problem can be solved analytically only 
in some special configurations. For example, if |V$| <C |g e 
everywhere in the system, a full perturbative analytical so- 
lution can be found (Milgrom 1986; Zhao & Tian 2006), but 
this is physically relevant only for the internal potential of 
low-density dwarf spheroidal galaxies orbiting more massive 
galaxies such as the Milky Way. 

In the Milky Way, if |g e | < a < |V$| (e.g. near the 
Galactic center), p m 1, and the whole dynamics is purely 
Newtonian. If |g e | <C |V$| <C ao (e.g. at 50 kpc from the 
Galactic center), then the dynamics is in the MOND regime, 
but the external field can be neglected, leading to a po- 
tential given by Eq. |T} up to a constant. Finally, when 
| \7<3> | ^ |g e j <c ao in the outskirts of the system, p tends 



to its asymptotic value p(\g e \/ao) — |ge|/ao, and the grav- 
itational potential is thus Newtonian with a renormalized 
gravitational constant. Then, contrary to the potential of 
Eq. Q, it is possible to escape from such a potential. This 
saturation of the /i-function is the essential piece needed to 
get a finite but nevertheless boosted escape speed (compared 
to its Newtonian value). 

This can be understood intuitively as follows. At the 
phenomenological level, the MOND force deriving from 
Eq. (JTJ can equivalently be interpreted in Newtonian gravity 
as providing an universal profile of dark halos. Explicitly, the 
effective dark mass M e fr at radius r 3> (GM/ao) 1 ^ 2 reads 
M off (r) ~ M X [(a Q ) 1/2 r/(GM) 1/2 - 1] and thus diverges 
with r. Interpreted in Newtonian gravity, MOND therefore 
predicts that any isolated baryonic mass M creates an ef- 
fective dark halo of infinite mass with density profile of the 
form p c g oc M e ff(r)/r 3 . The external field effect however 
saturates the /^-function at large distance, thereby leading 
to a potential which is logarithmic only on a finite range 
(GM/ao) 1 ' 2 -C r<r e , where r e is the transition radius at 
which the external field equals the local one: |g e | = |V$|. 
The external field thus truncates the effective dark halo at 
r ~ r e , yielding a finite effective dark mass surrounding the 
Galaxy, which may explain the oberved value of the local 
escape speed (Eq. [5}, depending on the actual value of g e . 



3.2 The escape speed 

To simplify the problem, we consider herafter the spherically 
symmetric case. The escape speed w osc at radius r is defined 
in a spherical potential as 

^LM = #(oo) - $ (r). (4) 

We assume that the spatial variation of the external field 
is so small compared to the size of the Milky Way that it 
is approximately constant over all space, so that $(oo) is 
well-defined. We use the assumptions and values from the 
best analytic MOND model of the Milky Way (Model III 
of FB05). This means that we assume that all the baryonic 
mass is confined to r < Ro , where Ro = 8 kpc is the galacto- 
centric radius of the Sun, and that the asymptotic circular 
velocity is v x = 165 kms -1 . The total baryonic mass M 
of the Milky Way is related to the latter by means of the 
Tully-Fisher relation GMao = f 4 ,. Finally the interpolat- 
ing p- function is chosen to be given by p(x) = x/(l + x). 
In addition to providing the best analytic fit to the termi- 
nal velocity curve of the Milky Way (FB05), this function 
has been shown to yield excellent fits of rotation curves in 
external galaxies (Famaey et al. 2007). 

We then consider as a toy-model a centripetal external 
field acting as a constant force directed towards the Galactic 
Center, in order to estimate the order of magnitude of the 
external field needed to produce the observed escape speed. 
We thus solve Eq. © for g = |V$| (we also note g e = |g e |). 
Using spherical symmetry and hence Gauss theorem, this 
means solving: 

flQff (ff + ge) _ U^o /gv 

a + g + g e r 2 ' 

Then we compute the local escape speed 
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Figure 1. Full line: the MOND potential of the Milky Way for the 
parameters of FB05 i?o = 8 kpc and v^o = 165 kms" 1 , and an 
external field directed towards the Galactic center g e = 0.01 X ao- 
The potential is such that <E>(i?o) = 0, and it asymptotes to 
<3?(oo) = 1/2d^ sc (_Ro) with v esc (Ro) = 544 kms -1 . Long-dashed 
line: the corresponding MOND potential without the external 
field. Short-dashed line: the Keplerian potential with a renormal- 
ized gravitational constant to which the MOND potential asymp- 
totes. 



«esc(-Ro) = ( 2 / g(r)dr 



1/2 



(6) 



and find that in order to have a local escape speed in the 
range allowed by Eq. ([2]), with a median of 544 kms" 1 , the 
external field must be 



g e = 0.010t o ;E£l x a . 



(7) 



Note that, in addition to the intrinsic non-sphericity of the 
inner Milky Way itself, the real external field points in a 
definite direction, meaning that the escape speed contours 
are in reality non-spherical, depending on the angle between 
the internal and external force vectors. As a first order test 
of robustness, note that replacing |g+ge| = (ff + ffe) by, e.g., 
(g 2 + gt) 1/2 in Eq. J5j) only leads to a 20 kms" 1 difference 
in the computed local escape speed for an external field g e 
of order ao/100. 

The Milky Way potential (whose zero is fixed at Ro) ob- 
tained from Eq. <(Sj is plotted on Figure 1. Overplotted are 
the corresponding MOND potential without external field 
(asymptoting to Eq. [IJ, and the Keplerian potential with 
renormalized gravitational constant Gao/g e - The actual po- 
tential interpolates smoothly between these two behaviours, 
the transition occurring when the internal and external ac- 
celeration are of the same order of magnitude. 



4 ESTIMATING THE EXTERNAL FIELD 

In this section we focus on the various known sources that 
may produce the needed external field (Eq. [7J) in a MOND 
Universe. One must however keep in mind that when there 
is more than one attractor the MOND gravity is not simply 
the sum of the contributions. 

There are two obvious possible contributions: the An- 
dromeda galaxy M31, and the Large Scale Structure (LSS). 
Andromeda (e.g. Evans et al. 2000) is located at a distance 
d,A ~ 800 kpc and has a baryonic mass Ma ~ 10 11 Mq. This 
means that the gravitational force per unit mass exerted by 
M31 on the Milky Way is of the order of 



{GM A a ) l/2 



(1.: 



ao/100, 



(8) 



with some variations depending on the actual location con- 
sidered in the Milky Way. However, since both galaxies have 
a roughly similar baryonic mass, the MOND gravitational 
field created by M31 is always smaller than the one created 
by the Milky Way in the direction opposite to M31, mean- 
ing that M31 cannot be the only source of the approximately 
constant external field we are looking for. 

On the other hand, estimating the total contribution 
from all the relevant distant sources in the LSS combined 
with the one of M31 is an extremely arduous task (the 
MOND gravity not being the sum of the individual con- 
tributions), and would require a detailed knowledge of the 
precise geometry and mass content of the MOND Universe 
at large. However, one can estimate that a gravitational field 
of order ao/100 is roughly the MOND acceleration 30 Mpc 
away from a galaxy cluster with an asymptotic, isotropic, 
line-of-sight velocity dispersion of 500 kms -1 (or 130 Mpc 
away from a cluster with a 1000 kms" 1 dispersion). This 
means that the Virgo and Coma clusters could both individ- 
ually contribute external fields of this magnitude (Milgrom 
2002). 

Moreover, it has been known for twenty years (Lynden- 
Bell et al. 1988) that the Local Group belongs to a 
600 kms" 1 outflow w.r.t. the reference frame of the Cos- 
mic Microwave Background. This is presumably due to the 
local gravitational attraction of LSS, and mainly from the 
so-called Great Attractor region at (I ~ 325°, b ~ —7°) (see 
Radburn-Smith et al. 2006). Once again, precisely calculat- 
ing the value of this gravitational attraction at the location 
of the Milky Way in a MOND Universe is virtually impos- 
sible, but one can estimate its order of magnitude from the 
acceleration endured by the Local Group during a Hubble 
time in order to attain a peculiar velocity of 600 kms -1 : 

H x eOOkms" 1 ~ ao/100, (9) 

for Ho = 70 kms" 1 Mpc" 1 . This estimation is independent 
from the assumption that we are living (or not) in a MOND 
Universe. In conclusion, M31, the Virgo and Coma clusters, 
as well as the LSS at large, all provide comparable gravita- 
tional fields that precisely lie in the allowed range (Eq. [TJl in 
order to produce the right local galactic escape speed. 



5 DISCUSSION 

Using the best analytic MOND fit for the terminal velocity 
curve of the Milky Way (Model III of FB05), and consider- 
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ing a toy external field acting as a constant force directed 
towards the Galactic Center, we have estimated the external 
field in which the Milky Way must be embedded in a MOND 
Universe in order to produce the observed local escape speed 
(Sect. 2, Eq. [2}. We found that the needed external field is 
of the order of a o /100 (Sect. 3, Eq. 0. 

One might question how this result depends on the spe- 
cific choice for the interpolating function (i. As we explained 
at the end of Sect. 3.1, the escape speed is mostly determined 
by the value of the transition radius r e where the internal 
and external accelerations are of the same order of magni- 
tude. The fact that the needed external field is well below 
ao means that the internal field at r e is in the deep-MOND 
regime fi{x) = x for any realistic /i-function. However the 
/i-function plays a non-trivial role in determining the stellar 
mass-to-light ratio needed to fit the terminal velocity curve 
of the Milky Way. For instance, using the "standard" func- 
tion fi(x) = x/(l + x 2 ) 1 / 2 requires a higher mass-to-light 
ratio, leading to an asymptotic velocity Voc = 175 kms -1 
(Model I of FB05). We calculate that it typically leads to a 
30 kms -1 increase in the local escape speed for an external 
field of order ao/100. Such a variation is thus small com- 
pared to the error bars quoted in Eq. @, and this signals 
the robustness of our estimate w.r.t. the choice of /i. On the 
other hand, our Eq. ([5]) is only an approximation for a couple 
of reasons: it neglects the presence of the distant external 
source(s) (and hence the direction of the external field), as 
well as a term proportional to the gradient of /i in Eq. @- 
However, both effects are proportional to g e , meaning that 
our result does hold qualitatively as a first order estimate. 

We argue that ao/100 is a value roughly compatible 
with the gravitational attraction that Large Scale Struc- 
tures are expected to exert on the Milky Way (Sect. 4). We 
thus conclude that the local galactic escape speed can be 
explained within a modified gravity framework a la MOND 
without the need of huge quantities of dark mass on galactic 
scales. However, it is very plausible that taking into account 
the matter density from M31, the Virgo and Coma clus- 
ters as well as all the structures within the Great Attrac- 
tor in a MOND potential solver (Ciotti, Londrillo & Nipoti 
2006) could perhaps lead to a local gravitational field slightly 
above the range allowed by Eq. ([7J| (note that an important 
issue regarding this calculation might be the possibility of 
a scale-dependence of ao; see e.g. the values used in Nusser 
2002 and Skordis et al. 2006). Such an increase of the exter- 
nal field would slightly decrease the corresponding galactic 
escape speed, and might signal the need for small quantities 
of dark mass (in baryonic or non-baryonic form) at large 
radii from the Galactic center (e.g. in the form of 2 eV neu- 
trinos; see Sanders 2003; Angus et al. 2007). 

As a final note, we stress that there are good theo- 
retical reasons to impose the /i-function to lie within the 
range e < n(x) < 1 for all x, where e is a small dimension- 
less parameter. Indeed, if the /i-function asymptotes in the 
standard way at x — > (as n{x) ~ a;) then the energy of 
the gravitational field surrounding a massive body is infi- 
nite, as was already pointed out in Bekenstein & Milgrom 
(1984). Moreover, the most successful relativistic theory of 
the MOND paradigm to date (Bekenstein 2004), may exhibit 
dynamical singularities if the free function is chosen in such 
a way that the MOND phenomenology is recovered with an 
effective /i-function asymptoting to zero at x — > (Brune- 



ton 2007). These two singularities are manifestly avoided if 
one instead uses a slight deviation from Milgrom's first pro- 
posal, namely fi(x) ~ e + x (Sanders 1986; Famaey et al. 
2007). This e-term acts precisely as the toy external field 
that we have used in Sect. 3. Comparison with observed ro- 
tation curves, and especially from the last observed radius 
of NGC 3741 (Gentile et al. 2007) were used in Famaey et 
al. (2007) to constrain e to be lower than 1/10. The present 
analysis of the local escape speed actually enables to put 
a stronger bound on this parameter. If e were of the order 
of 1/10, it would saturate the /i-function too early and the 
predicted escape speed would be significantly lower than the 
one observed. As a conclusion, the value of the local galactic 
escape speed constrains e to be e $ 1/ 100. 
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